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A Recovery-Based A Posteriori Error Estimator 
for H (curl) Interface Problems * * 


Zhiqiang Cai ' Shuhao Cao * 


Abstract 

This paper introduces a new recovery-based a posteriori error estimator for 
the lowest order Nedelec finite element approximation to the iT(curl) interface 
problem. The error estimator is analyzed by establishing both the reliability 
and the efficiency bounds and is supported by numerical results. Under cer¬ 
tain assumptions, it is proved that the reliability and efficiency constants are 
independent of the jumps of the coefficients. 


1 Introduction 


Let (1 C K 3 be a bounded polyhedral domain with Lipschitz boundary. Let & = 
{OdA 1 t* e a partition of the domain O with each subdomain f lj being polyhedron. 
The collection of interfaces ( LT/li dQj)\dQ is denoted by 3. Assume that g and /3 
are piecewise, positive constants with respect to the partition &. We consider the 
following iT(curl) interface problem: 


Vx(/j 1 Vxu)+/9m = /, in O, 

uxn = g, on <90, 


( 1 . 1 ) 


where n is the unit outward vector normal to the boundary of O. This model 
problem originates from a stable marching scheme of the second-order hyperbolic 
partial differential equation on the electric field intensity u that is resulted from the 
Maxwell equations (e.g. see [19, 23]). The g is the magnetic permeability, and the 
/3 ~ -^2 + ^ is related to the dielectric constant e and conductivity a scaled by 
the time-marching step size At. The boundary data g is “admissible” in a sense 
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that we will elaborate when introducing the finite element approximation (2.15). 
Throughout this article, boldface letters stand for vector fields and spaces of vector 
fields, non-boldface letters stand for scalar functions and spaces of scalar functions. 

The variational formulation of problem (1.1) involves the Hilbert space iT(curl), 
which is the collection of all square integrable vector fields whose curl are also square 
integrable: 

H (curl; Q) := {v € L 2 (Q) : Vxv e L 2 {Q)}. (1.2) 

The right hand side data f depends on the original source current and on the electric 
field intensity at previous time steps in the time-marching scheme. In almost all 
relevant literatures, f is assumed to be divergence free. In this paper, we assume that 
f £ H (div), where H (div) is the analog of the iT(curl) space for the divergence 
operator: 

H (div; 11) := {v £ L 2 {Q) : V-v £ L 2 (Q)j. (1.3) 

For the finite element approximation to (1.1), Nedelec introduced the iT(curl)- 
conforming edge elements in [28], which preserves the continuity of the tangential 
components, and certain a priori error estimates can be established (e.g. see [27]). 
However, the electromagnetic fields have limited regularities at reentrant corners and 
material interfaces (see [15, 17]). Hence the assumptions for a priori error estimates 
fail, and this is where adaptive mesh refinement is introduced to perform local mesh 
refining process within the regions that have relatively large approximation errors. 

The a posteriori error estimation for the H (curl) problem in (1.1) with constant 
or continuous coefficients has been studied recently by several researchers. Several 
types of a posteriori error estimators have been introduced and analyzed. These 
include residual-based estimators and the corresponding convergence analysis (ex¬ 
plicit [3, 18, 14, 32, 11, 12, 39], and implicit [22]), equilibrated estimators [5], and 
recovery-based estimators [29]. It is interesting to note that there are four types of 
errors in the explicit residual-based estimator (see [3]). Two of them are standard, 
i.e., the element residual and the face jump associated with the original equation in 
(1.1). The other two are also the element residual and the face jump, but associated 
with the divergence of the original equation: V-(/3tt) = V /. 

The recovery-based estimator studied in [29] may be viewed as an extension of 
the popular Zienkiewicz-Zhu (ZZ) error estimator ([40]) for the Poisson equation to 
the H{ curl) problem with constant coefficients. More specifically, two quantities 
related to the solution u and V xu are recovered based on the current approximation 
of the solution in a richer recovery space. The recovery space in [29] is the continuous 
piecewise polynomial space, and the recovery procedure is done through averaging 
on vertex patches. The resulting ZZ estimator consisting of two terms is shown to 
be equivalent to the face jumps across the element faces. The element residuals are 
not included in the estimator in [29]. 

The purpose of this paper is to develop and analyze an efficient, reliable, and ro¬ 
bust recovery-based a posteriori error estimator for the finite element approximation 
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to the H( curl) interface problem, i.e., problem (1.1) with piecewise constant coef¬ 
ficients. Theoretically, the efficiency refers that the local error indicator is bounded 
above by the local error, the reliability refers that the global error is bounded above 
by the global estimator. The robustness refers that constants in the efficiency and 
the reliability bounds are independent of the jumps of the coefficients. 

The recovered-based estimator introduced in this paper may be viewed as an 
extension of our previous work in [7, 8] on the diffusion interface problem to the 
iT(curl) interface problem, which partially resolve the non-robustness of ZZ error 
estimator for interface problems. Specifically, we recover two quantities related to 
(a~ 1 \7xu and /3u in the respective iT(curl)- and iT(div)-conforming finite element 
spaces (the lowest order Nedelec and Brezzi-Douglas-Marini elements respectively). 
For discussions on which quantities to be recovered and in what finite element spaces, 
see [7, 8]. The resulting estimator measures the face jumps of the tangential compo¬ 
nents and the normal component of the numerical approximations to /U -1 V xu and 
(3u, respectively. Our study indicates that the element residual is no longer higher 
order than the rest terms in the estimator, the contrary of which is proved to be the 
case in the diffusion problem ([9]). As a result, the element residual using recovered 
quantities is part of our proposed error estimator as well. 

Theoretically proving a robust reliability bound for the iT(curl) interface prob¬ 
lem is much harder than that for the diffusion interface problem. This is because one 
needs to estimate the dual norm of the residual functional over the iT(curl) space. 
To overcome this difficulty, one needs to use a Helmholtz decomposition of the error 
(e.g. see [3]). For the iT(curl) interface problem, additional difficulty is that the 
decomposition has to be stable under a coefficient weighted norm. To obtain such 
a decomposition is a non-trivial matter, for the discrete level, a discrete weighted 
Helmholtz decomposition is studied together with its application and analysis for 
non-overlapping domain decomposition method in [25]. Here for the continuous 
version, under certain assumptions, we are able to accomplish this task through es¬ 
tablishing a weighted identity relating gradient, curl, and divergence of a piecewisely 
smooth vector field (see Lemma A.4), which is an extension to the technique used 
in [16]. The final decomposition result in our paper is similar to the one in [38], 
the coefficient distribution setting is slightly more general than the one used in [38] . 
Our proof uses a piecewise regularity result from [17] than directly building from 
extension in [38]. 

Another necessary tool for proving a robust reliability bound is a tweaked version 
of the Clement-type interpolation. We are able to extend naturally from the idea in 
[4, 31] for the vertex-based continuous Lagrange element to the edge-based Nedelec 
element in three dimensions. Moreover, our quasi-monotone assumption on the 
distribution of the coefficients is based on edges which is similar to that of [31] 
based on vertices in 2D, and our proof borrows the idea from [4]. This is the first 
Nedelec interpolation known to achieve such a robust bound. 

Moving onto the efficiency bound estimate, we prove the every part of the local 
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recovery-based error estimator can be bounded by the robustly weighted residual- 
based error estimator. For the local error indicator measuring the irrotational part 
(gradient part) of the error, using the weighted averaging technique, the constant in 
the bounds is independent of the coefficient jumps across the interface. For the local 
estimator measuring the jump of the numerical approximation to /i _1 V x u (weak 
solenoidal or curl part of the error), the degrees of freedom of the corresponding 
recovered quantity sits on the edge (lowest order Nedelec elements). Consequently, 
the averaging is performed within an edge patch in 3D, and this resembles the 
averaging of Zienkiewicz-Zhu (ZZ) error estimator on vertex patch in 2D. It is known 
that, if the averaging is done in vertex patches that span across the interface, ZZ 
error estimator (even correctly weighted) is not robust with respect to the ratio of 
the max/min of the coefficients on diffusion interface problem (e.g. see [7]). Here 
for the H( curl) interface problem, under the assumption of the quasi-monotone 
distribution of the coefficients again, we are able to prove that the ZZ type averaging, 
if carefully weighted, yields a robust efficiency bound with respect to the coefficient 
jump. This is a first known result as well. 

Numerically, we are able to show that the recovery-based estimator studied in 
this paper is more accurate than the residual-based estimator in [3] for several test 
problems. 

This paper is organized as follows. The variational formulation and the H (curl)- 
conforming finite element approximation are introduced in section 2. The explicit 
recovery procedures and the resulting a posteriori error estimators are discussed in 
section 3. In section 4, the reliability and efficiency bounds are proved along with 
the technical tools for analysis. The proofs of the bounds for the weighted Helmholtz 
decomposition is presented in the appendix if certain conditions are met. Finally, 
some numerical results for the benchmark testing problems are presented in section 
5. 

2 Preliminaries 

2.1 Notations 

Hereby we list some formal definitions concerning problem (1.1). The function space 
for the variational problem is: 

iT g (curl;fl) := {u € lT(curl;D) : uxn = g on 3D}, (2-1) 

equipped with the H{ curl) norm 

I Ml a (curl) = (|MIl 2 (Q) + ll^ 7x ' U llz, 2 (Q) ) • (2-2) 

The bilinear form of the variational problem to (1.1) is: 

A(u, v) := / (p~ 1 V xu ■ V xv + (3 u ■ v)dx, (2.3) 

J n 
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and the coefficient-weighted norm related to this problem is: 

|||'u||| 2 := A(u, u). (2-4) 

If a subscript is added for the weighted norm, it means the local weighted norm 
defined on an open subset O C fl: 

|||w|||q := / (^ -1 V xu -V xu + /3u ■ u)dx. (2-5) 

Jo 

In addition to the standard H{ div) space (1.3), we need the weighted version as 
well: 

.ff (diva; fl) = {v E L 2 (Q) : V-(au) E L 2 (Q) in II}. (2-6) 

Let Th = U{IL} be a triangulation of D using tetrahedra elements. The sets of 
all the vertices, edges, and faces of this triangulation are denoted by A £p, and 
F\ j, respectively. Denote the vertices, edges, and faces being subsets or elements of 
a geometric object M. by A//, (M), £h(A4), and Fh( Ad), where M. can be an element 
from the objects in the simplicial complex of the triangulation like a specified element 
K , or the whole boundary dQ, etc. For any vertex z E A4, let be the nodal basis 
function of continuous piecewise linear element associated with the vertex z. 

A fixed unit normal vector np is assigned to each face F. and a fixed unit 
tangential vector t e to each edge e. For any scalar- or vector-valued function v, 
define J't;J F = v~ — v + on an interior face F € Fh with a fixed unit normal vector 
rip, where = lim v(x + enp). Define {u} P = (u + + v~)/2 as the average on 

e->0± ^ 

this face F. If F is a boundary face, the function v is extended by zero outside the 
domain to compute and {u} F . 

The following algebraic identity is handy later in proving identities involving 
interfaces for any scalar- or vector-valued quantities a and b: 

IHIf = + HfWf- (2-7) 

Denote the diameter of an element K G Fh by hp and the diameter of a face 
F G Fh by hp. We assume that the triangulation Th is shape regular (see [13]), and 
this assumption holds for any tetrahedron during the local mesh refining process. 

The following notations serve as the languages to describe the local element or 
face patches. They will be used later in local weighted recovery procedure (Sec¬ 
tion 3), and in the proof of estimates for the weighted Clement-type interpolation 
(Section 4). 

For a face F G Fh, let ujp be the patch of the tetrahedra sharing this face F. 
Let cjk.f be the patch of the tetrahedra sharing a face with K. 

For an edge e E £/,, denote by 

Ue = U {KeTh:ee£h ( K ) } F 
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(a) For edge e with the unit tangential vec¬ 
tor t e , oj e = uf = i Ki, uj e , F = Ui =i F t . 



(b) ujf = K+ U K- with the 
unit normal vector tif- 


Figure 1: A dissection view of the local edge patch cu e , ui e ,F, and face patch ujp 


the collection of all elements having e as a common edge, where £h(K ) is the collec¬ 
tion of edges of the element K. For the edge patch u/ e , we define two //-weighted edge 
patches associated with an edge e 6 £h, which can be understood as the collection 
of the elements with the biggest/smallest // _1 on an edge patch, are referred to 

u> e = UA'eXe-^’ where l e = {K C u e : /%' = max 

and Q e = Ua- e x e ^D where Z e = {K C uo e : = ^inin 

Denote the union of the interior faces within an edge patch as follows: 

u e,F = Re¬ 

using Figure 1 as an illustration, = uf =1 Fi. Define the //-weighted 
interior faces as follows: 

u e ,F = U Fei e F ’ where l e = {F G F h (Q e ) : F C w e ,F}- 

Taking Figure 1 as an example again, if ui e = K\ U I\ 2, then ij e: p = F\ U F2 U F4. 

For an element K G Th, denote the patch of all elements sharing an edge with 
K by 

UK,e = Uk&i k e K ' where l K . e = {I< G % ■ K C u e with e G £ h (K)}. (2.11) 

Similarly, for a vertex z G Afh, denote by 

Wz = UjA-eT^: zeAT h (K)} K 


( 2 . 8 ) 

(2.9) 
patch of 


( 2 . 10 ) 
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the collection of all elements having z as a common vertex. For the vertex patch 
u z , the /3- weighted edge patch associated with an vertex z G A4 is referred to 

u z = Uh-cT K - where l z = {K C u z : /3 K = max f3 K }■ (2-12) 

A Clu z 

For an element K E 7^, denote the patch of all elements sharing a vertex with K 
by 

u K ,z = l )k&i K z K ' wliere Z K,z = {K G Th : K C Uz with 2 G AT h (K)}. (2.13) 

2.2 Finite Element Approximation 

The corresponding variational formulation of (1.1) is 

{ Find u G iT„(curl; fi) such that: 

(2.14) 

A{u, v ) = (/, u), VuG iT 0 (curl; fi), 

where (■, •) is the standard L 2 (n)-inner product. Because of fi and (3 being uniformly 
positive on the domain, the coefficient-weighted norm (2.4) is equivalent to the graph 
norm (2.2) for H( curl; fi). Moreover the bilinear form (2.3) is intrinsically coercive 
with respect to this norm. By the Lax-Milgram lemma, there exists a unique solution 
in Ug ( cur 1: fi) to problem (2.14) when the boundary data is “admissible”. 

The solution u of (2.14) is approximated in a iT(curl; fi)-conforming finite ele¬ 
ment space: the lowest order Nedelec finite element space MT>q (see [28]). On each 
element K, define 

AfT>o(K) = {p(x) G Pi{K) : p = a + bxx, a,b G M 3 }. 

The global finite element space NT>q is glued together through the continuity con¬ 
dition of H{ curl; fi): 

AfD 0 = {p G iT(curl; fi) : p(x)\ R G VDq(A) VA G T h }. 

For simplicity, we assume that the Dirichlet boundary data can be represented as 
the tangential trace of an AfT>o vector field, i.e., g = pxn on the boundary, where 
p G ATP 0 . The finite element approximation is 

{ Find Uh G AfT>o 0 H q ( curl; fi) such that: 

(2.15) 

A(u h , v h ) = (/, v h ), V v h G AfV o O iT 0 (curl; fi). 

The problem in (2.15) is well-posed in its own right. 

Before building the error estimator, we need an iT(div; fi)-conforming hnite el¬ 
ement space BDAii, which is the linear order Brezzi-Douglas-Marini face element 
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(see [6]). On each element K, define, in a way that leads to the local basis construc¬ 
tion, 


BDA / l 1 (K) = {p{x) £ Pi(K) : p = a + cx + Vx(bs), a,s £ M 3 , c £ M 1 , b £ B(K)}, 


where B(I \) the space of quadratic edge bubble functions in K. Similarly the global 
KDA'ti inherits the continuity condition from H (div; fi): 

BDMi = {pe H( div; 0) : p(x)\ K £ BDMi(K) VK £ T h }. 


Let K £ Th be an arbitrary tetrahedral element with vertices Zj, Zj, Zk, and zi, and 
let rii be the outer unit vector normal to the face Fi = con v(zjZkZi), opposite to 
the vertex z t . Let t t j be the unit vector of the edge orienting in the direction 
of Zj — Zi. The ATDq nodal basis function for the edge can be written as (e.g. 
see [33, 35]) : 

XiTlj XjTli 
Pea — , . 

fc'lj * TLj tij • Tlj 

where A n for n = i, j, k, l are the barycentric coordinates associated with the vertex 
z n satisfying A, + Xj + A& + A; = 1. The degree of freedom of A/"X>o can be then 
associated with each edge e £ £h, in that ip e satisfies 


(2.16) 


Fe ‘ tp |g/ 


1 

e’ 



ip e • t e / ds 


±<5 ee /, \/e'e£h(K), 


where 5 ee > is the Kronecker delta. The plus sign is taken when locally i e ’s direction 
coincides with the globally fixed t e d s. 

Similarly, we cook up a customized version of the 1SDM. i nodal basis function 
associated with the vertex Zj on face Fj as follows 


= 9 


^jtjj 

Tli * t, i 


- 3 


^k^ik 


- 3 


A iUi 


W'i ' tik 'R'i ' til 


(2.17) 


Now the degrees of freedom of KZAA4i can be defined using the first moment on Fj 
because '4 , f z .z :i satisfies: 



■n Fi )X Zm dS 


±5. 


jmi 


for m £ {j, k. /}, 


and ■ TtF' = 0 for any face F' £ Fh{K) other than Fj. Similarly, the plus sign 

is taken when locally the exterior unit normal n* to face Fj is in the same direction 
with the globally fixed unit normal n Fi . 







3 The Recovery-type Error Estimator 


There are two important physical quantities of interest: the magnetic field intensity 
and displacement current density which are related to the electric held intensity u. 
The magnetic held intensity at current time step is denoted by a, and displace¬ 
ment current density diluted by the time step size is denoted by r. For iT(curl) 
problem (1.1) that is time-independent, they can be represented by the following 

cr = xu and r = /3u, (3.1) 

then the partial differential equation in (1.1) can be rewritten as 

V x a + t = f. (3-2) 

By the assumption of the data f E H( div;P), it is straightforward to verify that 
the u, cr, and r lie in the following spaces 

u€H( curl;0), t E iT( div;fi), 
and a E iT(curl; fi) n iT(div/r; fi), 

where the weighted space iJ"(div//; fi) is dehned in (2.6). 

The relation (3.3) indicates the continuity conditions u, cr, and r must fulfill 
in the continuous level in (3.1) and (3.2). These continuity requirements not just 
come from the operator theory in Hilbert spaces, but also translate from the original 
Maxwell equations. For an arbitrary interface S within the domain, if there is no 
surface charge on S, it is well known from physics (e.g., see [27]) that 

[uxn] s = 0, [(rxn] g = 0, Ijjct ■ n]] s = 0, and [[r • = 0. (3.4) 

These zero-jump conditions are consistent with the continuity conditions for H (curl) 
and H (div), respectively. However, during rendering the continuous problem into 
the finite element approximation, numerical magnetic held intensity and numerical 
displacement current density, 

cr h '■= M _1 VXU), and r h := f3u h , 

violate the second and the last continuity conditions from (3.4), respectively. There¬ 
fore, two quantities are recovered in the respective H (curl)-conforming and H (div)- 
conforming finite element spaces using an explicit local weighted averaging tech¬ 
nique. Note that that the normal component of is a piecewise polynomial of 
degree one on each face. Consequently, we need to use either KZ3A4i or ~KT i for 
recovering displacement current density, instead of ~KT o. 
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3.1 Local Recovery Procedure 

We recover two quantities, a^ and t£, based on ah and Th through weighted aver¬ 
aging locally on edge and face patches, respectively. To this end, for a fixed interior 
face F C Fh, denote by K± the neighboring tetrahedra sharing this F as a common 
face. Recall that np is the fixed unit vector normal to the face F, let K + be the 
tetrahedron with np as its inward normal, and K- with np as its outward normal 
(see Figure 1). Let 


7 f 


- 1/2 

Pk+ 


- 1/2 - 1/2 
/%_ + PK+ 


and = 


3 1/2 

Pk+ 




(3.5) 


and jp = 1 — jp, Kp = 1 — Up. If F is a boundary face with its neighboring 
tetrahedron K_ , we set jp = k p = 1 and 7 f = k f = 0. 

The local averages a^F and t^f on face F are chosen using the weights above: 


^ h,P — ip'-’h 


respectively, where ahi<± = 07 I , and Th f ± {x) = lim Th(x + enp) on face 

F. The notation discrepancy in above construction is due to the fact, which is 
mentioned earlier in previous subsection, that 77 ’s normal component on each face 
is a linear polynomial, yet ah is a constant vector on a fixed K. 

Now we construct the recovered quantities a* h and from the above local av¬ 
erages of ah. and 77 as follows: 


a h x = 


> = £ £ 


dS 


e££h 1 FCuj, 


e,F 


&e,F | JF 


> cp e (x) 


and Th(x)= V V IrL [ (T h ,p ■n F )X z ds) ^ FtZ (x), 
FeF hZ&mF )V FUF ) 


(3.7) 


where G A/iDo and i^f,z £ 13DA4.1 are the nodal basis functions associated with 
the respective edge e and vertex z on face F (see (2.16) and (2.17)). 

By the construction of the basis functions in (2.16) and (2.17), we can see a^ G 
U"(curl; P) and G iT(div; P), respectively. The degrees of freedom of are the 
weighted averages of 77 := /3uh on a face patch. The degrees of freedom of <rj£ is the 
weighted averages of ah := /i -1 V xmj on selected interior faces in an edge patch. 

Now we may define the local error indicator r) 2 K = Vk j _ + Vk 0 + Vp p based on 
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these averages plus the recovery-type element residual: 


VK,± = 


Vk,o = 


li xl ' 2 (T* h -[i 1 / 2 Vxu h 

r l/ 2 r* h - P l ' 2 u h 


L 2 {K) 


2 \ 2 

v± = 2 ^ r)K ^ 

K&T h 


L 2 {K) 


Vo = 


J 2 Vk, o > 


KeT h 


(3.8) 


1/2 11 * 2 X—^ 2 

VK,R = Hk h K 11/ - pu h - Vxa h \\ LHK) , and rfa = ^ q KR . 

KeT h 

The global error estimator is defined by rj 2 = rj 2 ^ + rj^ + r/ R . 


4 Reliability and Efficiency Bounds 

This section studies the reliability and efficiency of the estimators defined in the 
previous section. The efficiency bound of the local indicator is established in section 
4.4. To prove the reliability bound of the global estimator, we need two tools: (1) a 
Helmholtz decomposition with weighted norm estimate section 4.1, for detailed proof 
under certain assumption please see Appendix A) that splits the error into two parts, 
and (2) a modified Clement-type interpolation (section 4.2). Under the assumption 
of a robust weighted Helmholtz decomposition exists, two quasi-monotonicity as¬ 
sumptions on the distribution of the coefficients, the reliability bound is obtained 
in section 4.3, and it is uniform with respect to the jumps of the coefficients. 

4.1 Helmholtz Decomposition 

For any vector in Ho (curl; 12), there exists an orthogonal decomposition with re¬ 
spect to the bilinear form A(-, •) (see [17, 20]). This weighted splitting was used in 
[3] for continuously differentiable p, and /3 to prove the reliability bound of a residual- 
based a posteriori error estimator. Here we first present the robust weighted splitting 
result as an assumption (Assumption 4.1), then in Appendix A we show the proof of 
a bound independent of the coefficient jump ratio, under certain assumptions about 
the geometries and the relations between coefficients. 

Define the X (12, (3) as the space of curl-integrable functions intersecting weighted 
div-integrable vector fields, and PH S (12, 0 s ) as the space of piecewisely continuous 
vector fields on each subdomain: 

X(Q, P) = H(div/3; 12) n H 0 (curl; 12), 

9 I (4.1) 

and PH S (Q,&) = {v£ L 2 (Q) : G H s (%), j = 1 

Assumption 4.1 (Weighted Helmholtz decomposition). We assume that for any 
v G Ho (curl ; 12) ; there exist if G Hq(12) and w G PH 1 (H, &) n X(fl,/8) such that 
the following decomposition holds 

v = w + Vif. (4-2) 
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Moreover, the following estimate holds: 


5^ p~ 1/2 Vw 

3 =1 


+ 

L 2 (flj) 



+ 

L 2 (£l) 


/3 1 / 2 V^ 


L 2 {Vt) 


< C\\\v 


(4.3) 


4.2 Clement-type Nedelec Interpolation 

Weighted Clement-type interpolation operators for nodal Lagrange elements are 
studied in [4, 31]. Stability and approximation properties of this type of operators 
are often used in proving a robust reliability bound for a posteriori error estimators. 
For Nedelec elements, the standard unweighted quasi-interpolations for Nedelec ele¬ 
ments are studied in [3, 14, 32], In [3], the author defines the edge degrees of freedom 
by averaging on a certain face where that edge lies, which is similar to the construc¬ 
tion of the Scott-Zhang interpolation operators. In [14], the averaging is performing 
on the edge patch consisting of two triangles in 2D. Following the idea of averaging 
on the weighted vertex patch in [4, 31], and extending the averaging technique on 
edge patch in [14] to the three dimensional case, we construct a weighted Clement- 
type Nedelec interpolation operator from Jig (curl ;fl) to A/Dq. If the vector field 
to be interpolated has regularity, then the approximation and stability 

properties of the interpolant are proved to be robust under the weighted norm, with 
the assumption that the coefficient is quasi-monotone in Assumption 4.4. 

First we define the standard Nedelec interpolation in any K G Th■ To make this 
interpolant well-defined and bounded, we need to restrict that the vector field to be 
interpolated on each element K lies in the space H l / 2+5 {K) for some 5 > 0, with 
its curl in L P (K) for some p > 2 (see [27] Lemma 5.38). 

Definition 4.2 (Nedelec interpolation). For any v\ R G H L ^ 2+S (K ) with Vxv\ R G 
L P (K), define the interpolation operator : H 1 ^ 2 ^ 5 (K) —> M'Dq on each element 
K G Th as follows: 


EL v I is 


E 

e£S h {K) 


a e (v) ipe, with a e (v ) = f 

1C J e 


V 


t P ds. 


Definition 4.3 (Weighted Clement-type Nedelec interpolation). For any v G iTo( cur l; ^)> 
such that v \ G H (curl ; ojk e )nPi/ 1 ( ujk ei &), define the weighted quasi-interpolation 

'^K,e ’ ’ 

operator U h -L 2 m^AfD 0 on each element K G Th as follows: 


II h v \ K = a e(v)cp e , with a e (v) 

e£S h (K) 



■ t e 


if e is an interior edge, i.e., the 1-dimensional Lebesgue measure measi(e n <9fl) = 0. 
If e G Sh{dM), then a e (v) = 0. 
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To establish the stability and approximation bounds for this interpolation uni¬ 
form with respect to // _1 , a quasi-monotonicity assumption is needed on the distri¬ 
bution of the coefficients associated with each edge patch Lv e in three dimensions, 
which is similar to those of [4, 31] associated with each vertex patch in two dimen¬ 
sions. The quasi-monotonicity, in layman’s terms, can be phrased as “for every 
element in an edge patch, there exists a simply-connected element path leading 
to the element where the coefficient achieves the maximum (or minimum) on this 
patch”. The following assumption is stated in a mathematically rigorous way to 
convey above idea. 

Assumption 4.4 (Quasi-monotonicity of the /a^ 1 in an edge patch). For each edge 
e E £h, if e is an interior edge, for every K C oj e , and every K' C u e , (i) assume that 
there exist a collection of elements uj' e = \dSfK, C uj e with Ki(K )e ) C ui e , such that 
Ki shares a face with /Q_i, and that nfc < /iff for all i = 1,..., l(I\, e), where 
Kq = K. If e is a boundary edge, for every K C u ; e \w e , and every K' C u e , (ii) 
assume that (i) holds, and the 2-dimensional Lebesgue measure meas2(<9u; e n<9fl) > 0. 

The assumption is phrase using tD e , the assumption remains the same if we switch 
uj e to cD e , and reverse the direction of the inequalities. 

If Assumption 4.4 is satisfied, the extended p-weighted patch for an element K 
is denoted as 

U K,e = K u e e£ h (K) U e UeS^A')^- (4.4) 

Remark 4.5. Assumption 4.4 (i) is weaker than the extension of the quasi-monotonicity 
assumption in [4], and is the equivalent to the extension of the quasi-monotonicity 
assumption in [31] from the vertex patch in two dimensions to the edge patch in 
three dimensions. Notice if Assumption 4.4 is met, then ui e is a simply connected 
Lipschitz polyhedron for any interior edge e. 
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(a) Quasi-monotone in the 
senses of the extension of [4] , 
the extension of [31], and As¬ 
sumption 4.4. 



(b) Quasi-monotone in the 
senses of the extension of [31] 
and Assumption 4.4, not in 
the extension of [4]. 



(c) Not quasi-monotone in 
any sense. 


Figure 2: Different scenarios of the coefficient distribution for /i _1 for an interior 
edge patch iw e , where the edge e is marked as red dotted vertical edge in each figure. 
The tetrahedra whose bases are marked using blue color in (a) and (b) consist the 
u e for this edge patch. 


The illustrations in Figure 2 show the difference and similarity between the 
Assumption 4.4 and the extension to those in [4, 31]. In Figure 2a, for any two 
tetrahedra in the edge patch u e , there always exists a monotone path connecting 
these made of tetrahedra, along which one tetrahedron shares one face with the next 
tetrahedron in this path. In Figure 2b, along the path from any tetrahedron in this 
patch to the one with the biggest coefficients /r -1 , the coefficients are monotone. In 
Figure 2c, the coefficient distribution of the checkerboard type is not quasi-monotone 
in any sense, and a Clement-type interpolation cannot achieve a robust bound (see 
[31, 36]), if the edge of interest is an interior edge of the triangulation. 

Theorem 4.6 (Approximation and stability properties). Under Assumption 4.4, 
the interpolation operator f][h in Definition 4.3, satisfies the following estimates: 




~ 1/2 (« - n, 


v 


L 2 {K) 


< C hf{ 


h 


and 


h 


^vxf.-n, 


L 2 (K) 


< c 


" 1/2 Vr .. 

L 2 (uJK,e) 

H~ l / 2 Vv 


(4.5) 


L 2 (uK,e) 


for all v E iTo(curl;Q) such that v\ E iT(curl; UK,e) i?), where 

the Jacobian Vn is defined piecewisely by Vn|^ := V(n|^.). 

Proof. To establish the inequalities in (4.5), let % and v% e be the averages of v 
over K and u e respectively, i.e., vk = |A"| -1 Ik v dx, an d = |w e | _1 f~ v dx for 
an interior edge e. Let v% e = 0 if e is a boundary edge. 

If e is an interior edge, we have the following standard approximation property 
(also known as Poincare inequality) thanks to the shape regularity of the triangula¬ 
tion T h , simply-connectedness of uj e for an interior edge e from Assumption 4.4, and 
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(a) Quasi-monotone in the sense 
of Assumption 4.4, the extension 
of [31], and the extension of [4]. 



(c) Quasi-monotone for the edge 
patch in the sense of Assumption 
4.4, not quasi-monotone for 
vertex patch for the black dotted 
vertex. 



(b) Quasi-monotone in the sense 
of Assumption 4.4 and the 
extension of [31], not in the 
sense of the extension of [4]. 



(d) Not quasi-monotone in any 
sense. 


Figure 3: Different scenarios of the coefficient distribution for for a boundary 
edge patch u e , where e is marked red, the boundary faces are marked yellow, and 
the coefficient in each tetrahedron is marked on its front faces towards the viewer. 


\\ v - vk\\ L 2 {K ) < Ch K \\Vv\\ L 2 {K) , and \\v - v^ e \\ L2(Uje) < Ch e ||'Vv|| L a (Se) ■ 

(4.6) 

If e is a boundary edge, the first inequality in (4.6) still holds. To get an equality 
similar to the second one, the fact that v £ iTo(curl; D) and £ PH 1 (oJK,e, &) 

implies v t e L £ H 1 (u e ), and v t e I = 0. The following Friedrichs inequality 

I CDg I C/CiJq I 1 C/u L 

holds (even if u e is not simply-connected as in the case of Figure 3b) 


\ v ' M L*(Z e ) < Ch e ||V(V ' t e )\\ L 2^ e ) < Ch e \\^v\\ L ^u e ) 


(4.7) 


The starting point of the proof is to split the error we want to bound into parts. 
On any K £ 7ft, it follows from the triangle inequality that 


«-n h v 


L 2 {K) 


< ||l> — Vk I 


L 2 (K) 


+ 


V K - Uh V 


L*(K) 
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The first term can be estimated using (4.6) first inequality. For the second term, 
since W^vk = vk (with slightly abuse of notation we can extend vk to whole edge 
patch by letting it be its value on K ), we have the following partition on the element 
K _ 

VK - n h V = Uh(VK -v)= - VzJ • t e ip e . 

e&£ h {K) 

Now applying the triangle inequality, and using the fact that < C\K\ 1 ^ 2 

(see the construction of ip e in (2.16)) yield 


VK ~ U h V L2{R) < Y I \(VK-V2 e )-t e (p e 


'l 2 (k) 


e£S h {K) 


■ Y \(VK~VzJ-t e \ ■ \\Ve\\ L 2 {K) <C\K\ 1/2 Y \( V K -VxJ- t e \. 
e&£ h (K) e&£ h {K) 


(4.8) 


To establish the estimate for | [vk — v^ e ) ■ t e \ for each edge, we consider three 
cases. The first case is that when K c u e , using the triangle inequality, the estimates 
in (4.6) and (4.7) gives 


\K\ 1/2 \(vk ~ « 2 e ) ' te\ = \\{v K - vz e ) ■ t e \\ L2{K) 

<\\{VK-v)-t e \\ L2(K) + \\(V-Vz e ) -te\\ L 2 {K) 

< II VK - V\\ L2(K) + || (« - V Se ) ■ t e || L2( ~ e) < C~^~ 

A l K 


p~ 1/2 \/v 


L 2 {Ue) 


Here the term in front of the last inequality is treated as \\(v — v^ e ) ■ t e \\ L2 ^ e ^ < 
ll« - ««j£ 2 ( 5 j e) for an interior edge, and ||(v - ■ t e || L 2 (2e) = ||v ■ for a 

boundary edge. 

The second case is that when K (jL uj e , yet I\ is adjacent to u e , and we denote 
the face they share as dl\ n dui e = F. The fact that the tangential component of 
v along the edge e is continuous across the face F, and v £ PH 1 (uK,ej &) implies 
that v ■ t e £ H l {K Uw e U F) (e.g. see [27] Lemma 5.3). To establish the estimate, 
we need a standard trace inequality for p £ F[ l {K U uj e U F) (e.g. see [34] Lemma 
3.2): 

IIpIIl 2 ^) < C* j/i F / ||p||i 2 (^/) + hp ||Vp|| LW } , (4.9) 

where K' can be either the element of interest K, or the element K as a subset of 
u e which is adjacent to K. 

Now it follows from the triangle inequality and shape regularity of the triangu¬ 
lation that 


\K\ 1/2 \(v k -vz e )-t e \ = \\(v K - v%J ■ t e \\ L2{F) 

— Ch^ \\{vk ~ v ) ■ t e \\ L 2( F j + Chj( \\{y - vz e ) ■ t e \\ L2{F) . 
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The first term in the (4.10) can be estimated using (4.9) and then (4.6) 


h-x \\(vk v) ■ t e \\ L 2( F ) 

< Ch]l 2 {h ~ F l/2 \\(v K - v) • t e \\ L 2 (I<) + h}p 2 ||V(t> • t e )\\ L 2 {K) } 

<C{\\v K ~ v\\ L 2 {K) + h K ||Vv|| i2(x) } < Ch K \\Vv\\ l2(k) . 

1 /2 

For the second term hfi ||(v — v^ e ) ■ t e \\ L 2 ( F ) in (4.10), using the same argument 
yields a similar estimate, except passing the trace inequality from the face F to the 
element K C oj e this time: 

IK V — ^2 e ) • te\\i,2(F) — Ch-ft ||Vv||_ t 2(£De) • 

Combining the two inequalities obtained above gives the following estimate for any 
e E S h (K) thanks to < /x- 1 : 


\K\ 1/2 \(v K - V£ e ) ■ t e | < C 


llR 

-1/2 

l l K 




_1 / 2 v 


v 


L 2 (uj e ) 


(4.11) 


The third case is that when K <jLu e , nor does K share a face with u) e . By As¬ 
sumption 4.4 there is a simply connected patch consisting of K \,..., Ku Ke \_ 1 along 
which the /x _1 is monotone. Separating the term of interest by triangle inequality: 


\(v K - Vz e ) ■ t e | < \(vk - V Kl ) ■ t e \ 4-b (v K , 




- Vr 


then each of the above terms can be proved yielding the same form of estimate 
in (4.11) by the same argument. This result, together with the representation of 


V K - U h v 


l 2 {k) 


in (4.8), implies the first estimate in (4.5). 


For the second estimate in (4.5), using the inverse inequality, the triangle in¬ 
equality, and U h VK = vk again, we have that 


vxrv 


L 2 {K) 


Vxl]>-vir) 

fi h v ~ v 


L 2 (K) 


< Ch~ K l 


< C h^ 1 
-i 


n iX v ~ v k) 


L 2 {K) 


L 2 (K) 


+ Ch K ||v — vk\\l*{k) , 


which, together with the first estimate in (4.5) and (4.6), implies the second estimate. 
This completes the proof of the theorem. □ 


Assumption 4.7 (Quasi-monotonicity of the (3 in a vertex patch). For any vertex 
z E Mh, assume that the (3 satisfies the vertex patch quasi-monotonicity condition 
in [31]: if z is an interior vertex, for any K Cw z , and K' C ui z , (i) there exist a 
collection of elements oJ z = U'LTA C ui z with Ki(k, z ) C uj z , such that Ki shares 
a face with Ki-\ and that f3K i _ 1 < (3k, for all i = 1,..., l(K , z), where Kq = K and. 
If z is a vertex on the boundary, for every K C ui z \oj z , and every K' C ui z , (ii) 
assume that (i) holds, and the 2-dimensional Lebesgue measure meas2(c£; z n<9fl) > 0. 
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If Assumption 4.7 is satisfied, the extended /3-weighted patch for an element K 
is denoted as 

UK,z = A UzeM h (K) u 'z V)zeAf h (K)Uz, 

For the (3 which satisfies the vertex patch quasi-monotonicity in Assumption 4.7, 
the robust Clement-type interpolation for the linear Lagrange elements results are 
already established in [4, 31]. In the three dimensional setting, one reason to study 
the Clement-type interpolation is that the standard linear Lagrange nodal inter- 
polant may not be bounded. Unless extra regularity is assumed (e.g. the function 
to be interpolated is in i7 3//2+e (fl), see [27]), the degrees of freedom for the Lagrange 
nodal interpolant may not be well-defined because is not continuously em¬ 

bedded into the continuous function space. 

For any ?/> E 77q(Q), let V’h be the weighted Clement-type interpolant of (/> 
defined in [31] associated with the coefficient (3. Under Assumption 4.7, the V’/i has 
the following properties: 


and 


(3 1/2 (ip -'iphiz ))A 
/3 1/2 V- fa) 


L 2 (K) 


< Cl hj{ 


f3 1/2 Vip 


L 2 {ujk,z) 


L 2 (K) 


< C2 


f3 1/2 Vip 


(4.12) 




for any vertex 2 : E Afh(K). 



(a) Quasi-monotone in the sense of 
Assumption 4.4, not in Assumption 

4.7. 



(b) Not quasi-monotone in any 
sense. 


Figure 4: Different scenarios of the coefficient distribution for f3 for an interior vertex 
patch ui z . (a) /3 = 100 in the tetrahedra whose faces are marked blue, (3 = 1 in the 
rest tetrahedra in this patch. The coefficient distrubtion is quasi-monotone for all 
interior edges within this patch, (b) f3 = 100 in the four tetrahedra sharing the blue 
faces, (3 = 1 in the rest tetrahedra in this patch. The quasi-monotonicity is violated 
for the red edges. 


Remark 4.8. Assumption 4.4 does not exclude the case when D = Q 1 UU 2 , and 
Ui is a Lipschitz polyhedron touching the boundary at one vertex V\ only, with 


18 


















= 100, and = 1- Assumption 4.7 prohibits the existence of this 

scenario. In this scenario, a robust Clement-type interpolation cannot be achieved 
for nodal Lagrange elements. However, Assumption 4.4 allows this kind of domain, 
in which all the edges on the dLli connecting the vertex V\ is an interior edge of 
the triangulation. A robust Clement-type interpolation using Nedelec elements does 
exist in this scenario. Please refer to the illustration in Figure 3c. 

Remark 4.9. Assumption 4.4 which states quasi-monotonicity for the edge patch, 
is weaker than Assumption 4.7 for the vertex patch. The reason is that Assumption 
4.4 allows the checkerboard pattern for a vertex patch. However, this vertex patch 
checkerboard pattern is excluded in Assumption 4.7. Please refer to the illustration 
in Figure 4. In Figure 4a, the coefficient distribution satisfies Assumption 4.4 for 
any interior edges within this patch, yet Assumption 4.7 is not met. 

4.3 Reliability 

Under the assumption on the distributions of the coefficients and the Helmholtz de¬ 
composition which is stable under the weighted norm, we prove the global reliability 
for the local recovery error estimator p defined in (3.8). 

For any vertex z G Mh\Mh{dd), denote by 

Fui z := |~—| [ V-(/ — rf } )dx 

l W z| Juj z 

the average of V-(/ — rff) over the vertex patch u z . For z G J\fh{dLl), F,' Uz := 0. Let 

H = ( Y ^d] with r \K,d = P~K /2 hK ||V-(/-t£)|| L 2 W 

\KcTh ) 

and H = (j2 Y ||V-(/ - r* h ) - F^ . 

\zeM h K<zu z ) 

The contribution from interior nodes in H is a higher order term since V-(/ — r£) G 
L 2 (n), and so is the contribution from boundary nodes if V-(/ — t£) G IP(Ll) for 
some p > 2. (see [9]). 

Theorem 4.10 (Global Reliability of rj). Let u and be the solutions of (2.14) 
and (2.15), respectively. Under Assumption 4.1, 4.4, and 4.7 , there exists a constant 
C > 0 independent of the jumps of the coefficients such that 

III 14 — u h\\\ < C (77 + H). (4-13) 
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Proof. Denote the error and the residual by 

e = u-u h and R(v) = (/, v) - (/r _1 Vxw fe , Vxv) - (/3u h , v), 
respectively. It is easy to see that 

A(e, v ) = R(v), Vv G ffo(curl; O) and R(vh) = 0, G A/"X>o fl i?o(curl; Q). 

By Assumption 4.1, there exists a decomposition of the error e G iio(curl;fi) into 
the sum of if G Hq(£1) and w G Pi/ 1 (Q. iJ 21 ) n X(Q,j3) such that 

e = w + Vip and |||e||| 2 = R(e) = R(w) + R(Vip). 


To bound the curl-free part of the error, let iph be the weighted Clement-type 
interpolant of ip defined in [31] associated with the coefficient /3. It follows from 
the fact that RfViph) = 0, integration by parts, the Cauchy-Schwarz inequality, the 
approximation and stability of the interpolation (4.12), and (4.3) that 


P(VV0 = P(V(V> - 4> h )) = {f- t£, vty - W,.)) + K - P u h, V(V> - 4’h)) 
= -(V-(/ - T%),ip- ip h ) + (jl - /3u h ,V(ip - ip h )) 


s E 

KeT h 


m,d h 


-i 

K 


P 1/2 {4> - ip h ) 


L 2 (uj k ,z) 


+ VK, 0 


p 1 / 2 v( , ip - 4> h ) 


L 2 (ujk,z) 


< c ^2 (jlKA + r/K, o) 
K&T h 


(3 1/2 Vip 


L 2 (uk,z) 


<C(H + -no) 


p 1/2 V-ip 


L 2 (0) 


<C(H + r, 0 ) 


^ 2 < 


L 2 (f2) 


To bound R(w), let £ A/^Do D i/o(curl;D) with f| /( defined in 

Definition 4.3. Using the fact that R(wh) = 0 and integrating by parts give 


R(w) = R(w — Wh) 


= (f - /3u h , w-w h )~ (i tr * h , Vx(to- w h )) + (a* h - n 1 V x , V x (m - w h )) 

= (/ - Pu h -X7xal,w -w h ) + (<x ^ -|i _1 Vxu k ,Vx(tD -i«fe)). 

Now, by the Cauchy-Schwarz inequality, the second and third inequalities in (4.5) , 
and (4.3), we have 


R{w) < ^ (rjh'j{h K ] n 1/2 (w-w h ) 
K&T h ^ 


£ 2 W 


+ VK,X 


^ 1/2 Vx(»-u,,)| i](k) ) 


<C ^ (??if,i? + r?^_L) fi V 2 V 




L 2 (uj k ) 


< C(r] R + r )±) 


^Vi 


jj, 'Via 


i 2 (o) 


< C'(??fl + ??±) 


/i _1 / 2 V x e 


z, 2 (n) 
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Combining the above two inequalities, we have 

[||e||| 2 = R(e) <C(r] + H) |||e|||, 

which implies (4.13). This completes the proof of the theorem. □ 

In the remainder of this section, we assume that additionally V-/ E L p (Q) for 
some p > 2, then the H in (4.13) may be replaced by H which is a higher order 
term. 


Theorem 4.11 (Global Reliability of r/). Under Assumption A.l, 4.4, and 4.7, 
there exists a constant C > 0 independent of the jumps of the coefficients such that 

|||it — Uh \\| < C (j] + H^j . (4.14) 


Proof. In the proof of (4.13), if furthermore the following orthogonality condition 
is exploited on the vertex patch for the weighted Clement-type interpolant (e.g. see 
[7] Section 4) 

(l, (V* - ^hi.z))\ z ) ^ =0, Af h \Af h {dQ), 

together with the fact that F UJz = 0 and ifh{ z ) = 0 for z E A fh{dQ), it implies 

- T h)^-^h) K 

KeT h 

= ~ ~ F u z ,W~'iph{z))\ z ) K . 

z£j\fh KCu z 


Now, a similar argument as in the irrotational part proof of (4.13) gives 
{f-r* h y^-if h ))<CH 


l 2 (Vl) 


The rest of the proof for (4.14) is identical to that of (4.13). 


□ 


4.4 Efficiency 

Even though in [3], the coefficients are assumed to be continuous, the proof they 
used to prove the efficiency bound (Section 4 and 5 in [3]) carries over to piecewise 
constant coefficients. At the same time, their choice of weight yields a robust bound 
with no dependence on the coefficients. In this subsection, we prove the efficiency 
of the recovery-based estimator (3.8) by bounding the recovery-based local error 
estimator by the residual-based local error estimator. 
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Let fh be the standard L 2 -projection onto KDAii. It is proved in [3] that there 
exists a positive constant C such that: 


Ch F 

Ch F 

and 




L2 , F) < III 11 - u h\\\l F + Y VKh 2 K || f - fh\\ 2 L 2 {K) , 


KCuj f 


/3 F 1/2 l/3u h ■ nj 


Cn][ 2 h K \\f - /3u h - Vx(/j 1 V >< « fc )|| L 2 (A:) < III™ - u h\\\ K , 


L 2 {F) 


< 


- u h ) \ L2 ^+ Y Pk^k \\V-f\\l H K ), 


KCujp 


(4.15) 


where the coefficients fip 1 and (3 F on face F are given by the arithmetic averages of 
/r -1 and /3 

Vf 1 = (Ma_ + /%+)/ 2 > and + Pk+)/ 2 ’ 

respectively. Next we move on to prove the equivalence. 

Lemma 4.12 (Equivalence of f?A',o)* There exists a constant c > 0 independent of 
the jumps of the coefficients such that for any K E Th¬ 


em, o< Y ^f" P F l,2 W u h ■ njj 

FcdK 


L 2 (F) 


(4.16) 


Proof, it suffices to show that tjk ,o can be bounded by the summation of the residual- 
based estimator within the local face patch. 

For any interior element K, we first use a partition of unity to bound the esti¬ 
mator t)k ,o by the fact that NT>q(K) C BDM-i(K). The difference of the weighted 
average and t > t is 


(Th-Th)\ R = Y Y ( r^T Y Th ’ F - T h,K) -npXzdS] i/> Ftl 

FcdK\dfl z£Af h {F) ' 1 1 


E E 

FcdK\dfl zCATh(F) 


1 — K 


K 


|T| 


W u h ' n l F X * dS ^F, 


Recalling from (3.5) that on each face F of element K, k f = Px?/{P F 2 + Pk?)i 
where K' is the neighboring element sharing this fixed face F with K. Since 

II-MlRa) < C \F \ 2 an d W^PzWlz^k) — c\K \ 2 > 

and using the following coefficient weight relation (3.5) on each face F: 


(1 - k£)/ 3 k 1/2 = 


< 


V2 
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the local error indicator tjk ,o has the following bound: 


Vk,o 


< 


< 


r 1,2 r* h - p 1/2 u h 

E E 


L\K) 


1 


F<ZdK\dn z<zM h (F) mv F 
1/2 


1/2 


^/\r* h -r h ) 

W u h ■ ds 


E c ( jf ) ||^f 1 / 2 i^ u * - n ij 


FddK\dn 


L 2 (F) 


L 2 {K) 

W^F^Wl^K) 


For any element with a boundary face, thanks to the setting for problem (2.15), 
that the Dirichlet data can be exactly represented by an Af'Do vector field’s tan¬ 
gential trace, the degrees of freedom on any boundary face do not contribute to the 
approximation error in that element. This completes the proof of the lemma. □ 

Lemma 4.13 (Equivalence of rj Under Assumption 4.4, there exists a constant 
c > 0 independent of the jumps of the coefficients such that for any K e Th 


c r/K.A. < 


E E * 

e££ h (K) F<Zu b ,f 


1/2 


1/2 n 

Pf 


Km Vxujxn], 


L^F) 


(4.17) 


Proof. The proof of this lemma uses the setting in the edge patch’s illustration of 
Figure la. The edge patch u e consists of 4 tetrahedra, and the following proof gen¬ 
eralizes without essential changes to the case when there are more than 4 tetrahedra 
in oj e . 

Without loss of generality, the element of interest K is assumed to be K\ in 
Figure la. First performing the partition of unity for (Jh,K = M ^ xu ^| , which 
is a constant vector and can be represented by AfT>o(K) vector fields: 


VK,± = 




L*(K) 


V 2 E (°h _ a h,K) ‘ te Pt 

ee£ h (K) 


L 2 (K) 


* E 

e££h(K) 



&h,K ) * 


L 2 (K) 


< - a h,K) te 

ee£ h (K) 


Pe\\L 2 (K) ■ 


(4.18) 

By the fact that \\pe\\i,'z(K) — C\K\ 2 , the rest of the proof is to establish the 
equivalence, for every edge e, of | (erf — cr htK ) ■ t e | with the coefficient-weighted 
tangential jump term in the residual-based estimator. 

For the rest of the proof let us assume the edge of interest is e in Figure la. 
Before moving on to different coefficient distribution scenarios in this edge patch, 
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first by the local recovery (3.7) and the A/"X>o basis function construction (2.16), it 
is straightforward to check that 


t e = 




e,F 


E . 

FCU} e F 


(' cr h ,F ■ te) dS. 


The hrst case is when a) e = K = then tD e ,F = F\ U F 4 . Using the geometric 
relation that for any v ■ t e = n F x(vxn F ) ■ t e if t e lies on the planar surface F , and 
the definition of the weighted average 07 ^ in (3.6), yields 


[Ph h,I< ) ‘ t e 



iPh,Fi 


&h,K) ' te dS 



n Fi X (( (T h)Fi - cr h ,K ) x n Fi ) ■ t e dS 
(!-7 F i )K^~ 1 Vxu h) xn l Fi ’ ( t e xn Fi )dS. 


(4.19) 


By the coefficient weight dehned in (3.5), for F\ we have 


1 / 2 F , 

Vk t 1 “7 fJ = — 


K ~ 


1 


%/2 


1/2 


^ 1/2 + ^ /2 -(^+^) 1/2 
By Cauchy-Schwarz inequality and the triangle inequality 

Fk 2 \(*% ~ <*h,K) • t e | < ^ |Ti | 1/2 gffK^VxufcJxnlj 


2 w e i?i 
1 C) * 1 *G{1,4} 


LHFi) 

(4.20) 

Then using the shape regularity of the mesh, i.e. |i 7 j| 1 / 2 |_A| 1 // 2 |u5 e ^| _1 < C iijl 2 
for any i 7 ) in this edge patch, we have 


- a KK ) ■ t e 


L\K) 


< 


h F 2 1 Vxw h )xnJ J 


FCuJe } F 


L 2 (F) 


.(4-21) 

A variant of the hrst case is that K = K 1 C uj £t p. Assume u) ejF = K± U A' 2 , then 
Q e>F = F\ UF 2 U A 4 . By the dehnition of uj e , F in (2.9), = min,- ^...4 /j^.’ . 

The proof of the bound (4.20) for this variant shares almost the same argument 
with above, except there will be one extra term comparing to (4.19), and it can be 
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rewritten as follows: 


1 


\We,F\ JF 2 

i r 


\We,F\ J F 2 l 

1 


&h,Fo &h,K) ‘ t e dS 

(°7 i,F 2 ~ a h,I<2 ) + (O’/I,~ O’/i./f) 


• £,= d5 




+ 


" r [ (1-7‘Vxukjxn] ■ (t e xn F2 )dS 
,f\ Jf 2 

—r/ I(ii _1 Vx«;,)xnL • (t e xn Fl )dS. 

e,F I JFo 1 


(4.22) 


Using the the shape regularity of the edge patch (c|^ 2 1 < |Ti| < C*||), and the 
fact that 


d'K (1 7F 2 2 ) — Ma 2 (1 7f 2 2 ) < Hf 2 ’ an d d'K ~ _i/9 _i/2 — 2 / J 'i4 

Pk +Vk 2 


1/2 


1/2 


we reach the following estimate 


<- 


1/2 

/% 

C 


-7 [ (&h,F 2 ~ <Th,K) ■ te dS 

\^e,F\ JF 2 


\ Fi \ 1/2 til? Ka 6 1 Vxu h )xn]] F . 


l w e,F| . 


* 6 ( 1 , 2 } 


(4.23) 


i 2 (U) 


Thus the estimate (4.21) follows. If cD e contains more elements, the same ar¬ 
gument with above applies, with all the unweighted extra terms involve only the 
interior faces of Q e . This completes the proof for the first case. 

The second case when K = K\ (jL £) e , yet K is adjacent to Q e . Assume K 2 = u) e , 
i.e., u e ,F = F\ U F' 2 - A similar split as (4.19) applies 


(07 &h,K) ' te 



&h,K ) ' t e dS. 


The F\ term can be estimated the same with (4.20). The F 2 term can be rewritten 
as (4.22). This time we use > fij^ = minj = i r .. ) 4 /x^, this implies 


1/2 


1/2 , 

i l K + l J K 2 


-1/2 < 2 df. 


thus the estimate (4.23) follows, which, under some backtracking, confirms the va¬ 
lidities of estimates (4.20) and (4.21). If u) e contains more elements than I\ 2 , same 
argument applies as long as /i^- 1 > /j 7 ] and the shape regularity holds for the edge 
patch of interest. This completes the proof for the second case. 
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The third case is that K = K\ <£_ u; e , nor is K neighboring to Q e . Assuming 
u} e = K 3 , then uj e) p = F 2 U T 3 . The same split with (4.19) applies, but this time on 
face F 2 and F 3 , 


(&h ^h,K ) ‘ te 



(Vhfi 


&h,K ) ' te dS. 


In this lemma, Assumption 4.4 holds. Without loss of generality, we assume the 
monotone path from K = K\ to K 3 is through K 2 . The F 2 term can be estimated 
exactly like previous case, because — /%!. = m i n i=i ...,4 For the -^3 

term, using the same trick as (4.22) yields: 



07 x,k) 


t e dS 


By the quasi-monotonicity of the coefficient on this edge patch again, we have 

A/ 2 fl ^ „l/ 2 n ^ „l /2 n „ A ,. 1/2 - 2 

Fa 


/ 4 H 1 - 'Yft) < /^}/ 3 2 (l - 7^ 3 3 ) < and < — 


Ts's> 

3 


4/2 

P 3 


i/2 . 

+ ^K 3 


— 1/2 < . 


therefore, the estimate for the T 3 term is similar to (4.23), with one extra face in¬ 
cluded due to the fact that the inequality is passed through an intermediate element 
along the monotone path 


<- 


1/2 

l J K 
C 


rw—7 f {cTh, f 3 ~ &h,K ) • t e dS 

\^e,F\ JF 3 


\ Fi 1 1/2 VnKl* X V xuh)xn} Fi 


16 ( 1 , 2 , 3 } 


(4.24) 


L 2 (Fi) 


Consequently, estimates (4.20) and (4.21) follow for the third case. If the reader 
walks through the proof, one will find that more tetrahedra being contained in Q e 
than 1 does not change the essential part of the proof because of the existence of 
the monotone path. This completes the proof of the lemma. □ 


Theorem 4.14 (Local Efficiency of tjk)- Under Assumption 4-4> there exists a 
constant c > 0 independent of the jumps of the coefficients such that for any K E Tf: 


crj K < |||u - u h \\\ UKF +osc (f,n,P;u K , F ), 


(4.25) 
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where osc(/,/r, /?; ojr,f) is the oscillation of the data within ujk,f 


OSC (f,H,p-,U Kl F) = < Y {pK h2 K IIV • f\\ 2 L 2 {K ) + Mch 2 K \\f ~ fh\\ 2 L 2(K)) 


1/2 


K(Zlok f 


Proof. By the residual-based estimator local efficiency estimate (4.15), Lemma 4.12 
and 4.13 which show the recovery-based ijk ,o and rjx ._l can be bounded the face 
jumps in the residual-based estimator, it suffices to show that the local recovery- 
based residual term is locally efficient. Applying the triangle inequality for tjk,r 
gives: 


'IK,, 


which, together with a standard inverse inequality and (4.15), shows that 


cr]K,R < |||m - u h \\\ K + rj K)± . 


This completes the proof of the theorem. 


□ 


5 Numerical Experiments 

This section reports numerical results of our estimator on several three dimensional 
iT(curl) interface test problems. 

The numerical tests are implemented under ZFEM (see [10]) framework in MAT- 
LAB. Initial meshes are generated by the MATLAB built-in DelaunayTri and 
distmesh (see [30]). At each iteration, let Sh be a subset of Th whose elements 
satisfy 

Y 7 1k- 9 Y 

KeS h KeT h 

where the rjx is evaluated using the recovered quantities computed by weighted L 2 - 
projections through multigrid 17(3, 2)-cycle iterations. This procedure is analyzed 
in [37] for diffusion problem, and is proved to be equivalent to the local weighted 
averaging. The marking parameter 6 is chosen to be 0.2. All elements in S} L are 
refined locally by bisecting the longest edge, and some neighboring elements of Sh 
are refined to preserve conformity of the triangulation. 

To measure the global reliability of the a posteriori error estimator, we show 
comparisons of different measures in the each example’s table of comparison, n is 
the number of levels of refinement. The N n the dimension of £h, n in the n-th level 
triangulation, in our case, it is the number of degrees of freedom. The effectivity 
index for each estimator at the n-th level is: 

eff-index := ^ 

|||^ — ’ u '/i,n||| 
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where r] n is the error estimator, and Uh n is the finite element approximation at the 
n-th level of triangulation. 

The orders of convergence are computed for both r/ and |||w — Uh |||. r v and r err 
are defined as the slope for the line of rj n and |||u — •u/ l>n ||[ in the log-log scale plot, 
such that 


hi% ~ -r n lniV n + ci, and ln|||u - u h , n ||| ~ -r ei r lniV n + c 2 . 

In the convergence rate plot, the log of degrees of freedom is the horizontal axis, 
and the log of the error/estimator is the vertical axis. The order of convergence is 
optimal when r v and r eiT are approximately 1/3. 

In first two examples with known true solutions, the adaptive mesh refinement 
procedure is terminated when the true relative error 

rel-error := |||w — Ufc|||/|||w||| < Tol. 

For comparison, numerical results involve some of the following error estimators 
other than the recovery estimator in (3.8): 

1. The residual estimator in [3]: 

VK,Res = HKh 2 K \ \ f ~ /3u h - V X (/i -1 Vx« h ) ||^ 2(x) + Pj?h 2 K ||V-(/3u h - f)\\ 2 L 2( K ) 

+ + t J ‘ F ||[(M _1 VXU h )xn]] F ||^ 2(F) ) ’ 

F&F h (K) 

(5T) 

and rff ies = YlxeT h/v R es i where p F l and (3f are the arithmetic averages of /i 1 
and /3, respectively, on elements sharing the face F. Note that this estimator 
is weighted appropriately and may be viewed as the extension of the residual 
estimator in [4, 31] for the diffusion interface problem to the iT(curl) interface 
problem. 


2. The Zienkiewicz-Zhu (ZZ) based error estimator in [29] using the coefficient- 
weighted norm (2.4): 


Vk.zz ~ 


p 1//2 7£_l(Vx Uh) — p ^Vxuf, 


2 

L\K) 


+ 


/3 1/2 n 0 (u h ) - p 1/2 u h 


(5.2) 


L 2 (K) 


and rj' zz = Vk zz • Both recovered quantities 7 H_l(V xu/,) and TZo(uh) 

are in the continuous piecewise linear vector fields space Vi := {p E H l (Q) : 
p\ K E P\(K)}, and their nodal values at any vertex z E A// are: 


n±('Vxu h )\ z 


w 2 


V x Uh dx and 7Zo(uh) I 



Uh dx. 
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3. The Zienkiewicz-Zhu (ZZ) flux based error estimator in [29] with weight suited 
to the coefficient-weighted norm (2.4): 


VK,ZZ,f 


/x 1 ^ 2 TZ_i(n 1 Vxit/j) — [i 1 / 2 S7xui 1 


L 2 {K) 


+ 


p- 1/2 n 0 (pu h ) - p 1/2 u h 


L 2 (K) 


(5.3) 


and r] 2 zz j = YlxeTh^K zz f Both recovered quantities 77 j_(/r 'Vxjij,) and 
TZo(/3uh) are in Vi as well, and their nodal values at any vertex z 6 Mh are: 


TZ±(n 1 Vxu/j)| = 


w 2 


(x 1 Vxm/j dx and TZo(/3uh)\ = 


COz 


/3iih dx. 


Example 1: This example is adapted from a benchmark test problem (see [7, 
8, 26]) for elliptic interface problems. The computational domain is a narrow slit 
along z-direction: Q = (—1, l) 2 x (—5, S) with 6 = 0.2. The true solution u is given 
in cylindrical coordinates ( r,6,z ): 

u = Vi/> = V(r a <l>(0)), 

where (j>(6) takes different values within four different subdomains while being glued 
together using continuity conditions that is firstly invented in [26]. The fx = 1, and 
the (3 is given by 


(3 


R in (0, l) 2 x (—6, 5) U (1,0) 2 x (—5, <5), 

1 in Q\ ((0, l) 2 x (-5,5) U (1,0) 2 x (-5,5)). 


Here we set parameters a, R to be 


a = 0.5, R « 5.8284271247461907. 


In this example, the tolerance is set to be Tol = 0.1. The numerical results of 
example 1 are in Table 1. It shows that to achieve approximately the same level of 
relative error, the number of degrees of freedom needed in the mesh refined by the 
local indicator r]K,zz or rj k,zzj requires more than twice than the other two. 

The adaptively refined mesh generated by each estimator can be found in Fig¬ 
ure 5. The tendency of r\zz or rjzzj to over-refine those four interfaces is due to the 
fact that recovered quantities enforce unnecessary extra continuity conditions of the 
true quantities. For example, xuf t ) and 77o(/3ith) in (5.3) are in i/ 1 (H), 

yet for the true solution u, /i -1 V xu £ H( curl; fi) and (3u £ H( div; fi) in (3.4). 

Overall, the recovery-based error estimator and residual-based error estimator 
lead to the correctly refined mesh, and the recovery-based one performs more con¬ 
vincingly showing a less oscillatory convergence, achieving the same level of relative 
error in fewer iterations. More importantly, it exhibits a better effectivity index. 
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(a) Refined Mesh based on rjzz,K■ (b) Refined Mesh based on rjzz,f,K- 




Figure 5: Mesh result of Example 1 


Table 1: Comparison of the estimators in Example 1 



n 

# DoF 

rel-error 

eff-index 

r v 

^"err 

Vzz 

30 

39391 

0.0898 

1.309 

0.128 

0.235 

7 lzz,f 

26 

44111 

0.0837 

1.495 

0.243 

0.253 

VRes 

24 

18832 

0.0873 

1.749 

0.257 

0.301 

V 

18 

18649 

0.0886 

0.820 

0.299 

0.303 


Example 2: This example is in the numerical experiments section of [24], The 
domain is = B2 = {(x, y, z) : x 2 + y 2 + z 2 < 2 }, and the coefficients are given by 

(fi = m = 1, /3 = 1 in Si = {(x, y, z) : x 2 + y 2 + z 2 < 1}, 

|/i = /j /2 = 10 6 , j3 = 1 in Q\B\. 

The true solution u is given by yu\ in B 1 , and jiu-i in £l\B 1 . For the explicit 
expression please refer to [24]. The Tol = 0.2 in this example. 

In this example, the element residual term r)R in (3.8) is not a higher order term 
(see Figure 7a). The red dashed line is a reference line of a constant multiple of 
(^DoF) -1 / 3 . The numerical results of example 2 are in Table 2. The adaptively 
refined mesh of each estimator can be found in Figure 6 . 

The refined meshes based on rjK,Res, and rjx respectively are visually similar, the 
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r]K,zz and t]k,zzj tend to over-refine the region where the local coefficient-weighted 
error is not significant yet xu is discontinuous across the interface. 


Table 2: Comparison of the estimators in Example 2 



n 

# DoF 

rel-error 

eff-index 

r v 

^err 

r/zz 

18 

200692 

0.199 

4.077 

Not converging 

0.118 

Vzzj 

20 

99794 

0.193 

2.527 

0.839 

0.169 

VRes 

9 

63405 

0.186 

1.761 

0.273 

0.236 

V 

8 

52287 

0.193 

1.079 

0.282 

0.251 



(a) Relative error distribution of refined mesh (b) Relative error distribution of refined mesh 
based on rjK,zzj cut on y = 0, it can be ob- based on tjk cut on y = 0, the local errors are 
served that the local errors are not evenly dis-more evenly distributed than the mesh refined 
tributed. based on 


Figure 6: Relative error distributions result of Example 2 

Example 3: This example is a widely-used test problem examining the perfor¬ 
mance of adaptive mesh refinement procedure for Maxwell equations (e.g. see [23]). 
The true solution is unknown and not smooth. The homogeneous Dirichlet bound¬ 
ary condition is enforced, together with a constant source current f = (1,1,1). The 
coefficients are given by: 

\ H = 1, P = 1 in D c , 

= 1, /3 = 100 in fi\n c 

where fl c is {(x,y, z) : |x|, |y|, \z\ < ^}. In this example, we cannot compute the 
true error, hence we set the stop criterion to be r/ < Tol with Tol = 0.16 at the 
n-th level of triangulation. 

This example illustrates two important aspects: (1) the element residual is in¬ 
dispensable in the error estimator in the pre-asymptotic region; (2) the iterative 
refining procedure using the residual-based estimator engages much more degrees of 
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(a) Example 2: convergence of 77 , comparing to (b) Example 3: convergence of rj versus the 
the element residual term r/R. recovery-based only rm ec = (r]± + ^o) 1 • 

Figure 7: Convergence Results of Example 2 and Example 3 
Table 3: Comparison of the estimators in Example 3 



n 

#(DoF) 

Estimator 

r v 

VRes 

21 

219993 

0.156 

0.328 

V 

14 

61302 

0.152 

0.346 

iRec 

8 

15672 

0.159 

0.230 


freedom than the one using the recovery-based estimator, when same stop criterion 
is used for both. 

The recovery-based error estimator shows an optimal order of convergence, which 
is ~ (^DoF) -1 / 3 , i.e., r, ? ~ 1/3, if the local error indicator includes the element 
residual rjx '■= ( Vk± + Vk 0 + Vk r) 1 ^- ^ the element residual is discarded, i.e., the 

pure recovery-based estimator i]K,Rec ■= (Vk,± + Vk, o) V2 is used as the local error 

/ ’ \ 1/2 

indicator, the order of convergence for r]R ec := i Vk Reel no ^ optimal (see 

y K&T h ’ J 

Figure 7b, and Table 3). 

From the first two examples, we learn that the effectivity index of the recovery- 
based estimator is in general two times as effective as that of the residual-based 
estimator. For problem with an unknown solution which is quite common originated 
from some real world applications, when setting the stopping criterion using the 
global error estimator, the number of degrees of freedom using the residual-based 
error estimator is (eff-indexRes/eff-indexRec) 3 as much as that using the recovery- 
based error estimator (see Table 3). 
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A Weighted Helmholtz Decomposition 


Here we establish a weighted Helmholtz decomposition in light of [16, 17] tailored for 
the H{ curl) interface problem. The following assumption is needed to guarantee 
that such a decomposition exists with the constant in estimate (4.13) is indepen¬ 
dent of the jumps of the coefficients. In other words, the constant in the estimate 
depends on the jump size of the product of two coefficients, and the geometries of 
the interfaces as well. 

Assumption A.l. (i) The domain H is assumed to be convex, simply-connected, 
and that no three or more subdomains share one edge from the triangulation o/H. 
(ii) The coefficients /j and /3 are assumed to satisfy: C m in < pjftj < C max , where 
C m in and C max are two constants independent of the jumps of p,j and (3j, or flj on 
each subdomain Qj. 

Firstly, we define some additional function spaces, along with the X(Q,(3) in 
(4.1), relevant to the weighted Helmholtz decomposition as follows: for any piecewise 
constant a = a.j in Qj : 

iT 0 (curl; SI) = {u £ iTo(curl; P) : Vxu = 0 in fi}, 

PC°°(H, &) = {v £ L 2 (H) : t,| n . € C 00 ^), j = 1,..., m}, (A.l) 

and PH S (Q, 3?) = {«£ L 2 (Q) : u| Q £ j = 1, ..., m}. 

O 

It is well known (see [21]) that the kernel of curl operator, H o(curl;H), is 
characterized by the gradient held in a simply-connected domain: 

O 

Lemma A.2. If Q is simply-connected, for any u £ Lfo( cur b^)) there exists a 
unique function if £ H ( \ (Q) such that u = Vif. 

O 

Since A(v,Vif) = {j3v, 'S/if), the orthogonal complement of i3"o(curl; P) = 
V77 q(H) with respect to A(-, •) is 

{v £ ffo(curl;fi) : (Pv,Vip) = 0 \hf £ H ( ](Q)} c X(Q,p). 

To construct a weighted Helmholtz decomposition tailored for the interface prob¬ 
lem, an analysis of the structure of X(£l,f3) is necessary. Before tackling this, the 
following lemma from [17] is needed: 

Lemma A.3. X(£l,a) n PC 00 (Cl, &) is dense in X(P,a) n PH 1 (Cl, CP) in the 
following norm: 

ll^llx(n) = ll v lli, 2 (n) + ll Vx,y lli 2 (n) + ll v '( m; )llL 2 (n) ■ 

Now we move on to prove the norm equivalence for certain piecewise H 1 -ve ctor 
fields using the density argument of Lemma A.3. 
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Lemma A.4 (Norm equivalence for piecewise smooth vector fields). For all v E 
X(Q, a) n PH l (Q, 0 21 ), the following identity holds: 

m r. 

E/ a\Vv\ 2 = / (a|Vxu | 2 + cT 1 ] V-(cru)| 2 ) . (A.2) 

j J i2 j J S2 

Proof. By Lemma A.3, it suffices to establish identity (A.2) for any 0 E A(fl, a) n 
PC°°(Q, 0 s ). To this end, using a local identity —A <f> = Vx(Vx0) — V(V-0) and 
integrating by parts on each subdomain Qj twice give: 

m r. m „ 

E/ «|V0 | 2 = E/ a j (|Vx 0| 2 + |V- 0 | 2 ) + B 

j =i j =i 

m ~ 

with B = • V)0 + nx (Vx0) — (V-0)n'j • <f>dS. 

, i ■/an,- v > 

Now, it remains to prove that B = 0. On any polygonal face with normal vector 
n that is represented by the cartesian coordinates in the three dimensional space, 
rather than the local planar coordinates, 0 may be decomposed into the normal and 
tangential components as follows: 


0 = (0 • n)n + 0 t with 0 j = rax( 0 xn), 

which, in turn, implies 

(n • V)0 = ((n • V)(0 • n))n + (n ■ V)0t 
and V -0 = V - ((0 ■ n)n) + V-0y = (n • V)(0 • n) + V-0 t- 

By using the following identity (e.g. see [2]) 

V(a • b) = (a • V)b + (6 • V)a + ax (Vx b) + bx (Vxa) 

and noticing that n is a constant vector on a face, we have 

V(n • 0) = (n • V)0 + nx (Vx0), 

which yields the following by being projected onto each polygonal face 

Vj(n • 0 ) = (n • V)0t + nx (VX 0 ), 


(A.3) 


(A.4) 


(A.5) 


where Vj is defined as Vju = (Vn)j. It follows from (A.4), (A.5), (A.3), horno- 
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geneous boundary condition, and identity (2.7) that: 


B = Y'' / aj (Vt(w • 0) ~ (V-0t)?"i ) ' 0 dS 

j=i V y 

m r. 

= Y,j a i ( v T(n • 0) • 0T - (V-0 T )(0 • n)) dS 

j=l 3 9Clj 

= Y1 f (l v T(«0 • n ) ■ 0t1 f - I(V-0t)(Q! 0 • n)J F ) dS 

FC3 Jf 

= [ ( v t ([«0 • ™} F ) • 0 t + v - [ 0 t 1 f ( Q !0 • n ) 

fc3^ f 

+ V T (a0 • n) • [[0 T 1 F + V• 0 T [o:0 • n]] F ) dS 1 . 

Now B = 0 is a direct consequence of the continuity conditions for 0 E A(f2, a) n 

PC°°{n,&): 

[[0xnJ F = 0 and [[a0 • nJ F = 0 VF C 3. 

This completes the proof of the lemma. □ 


Remark A.5. Lemma A.4 is an extension to the Lemma 3.8 in [21] for Lipschitz 
polyhedron in the case when only homogeneous tangential boundary condition is sat¬ 
isfied for the vector field. It uses a similar argument to that of the Theorem 2.3 in 

[16] . In [16], no piecewise constant coefficients are involved, but the technique used 
shed light upon this kind of identity. The result in Lemma A.4 bears the same form 
with an identity valid for PH 2 {fit , IP) regular vector fields used in Lemma 2.2 in 

[17] . In the proof of Lemma A.4, we further exploit the density result in [17], which 
implies this identity in [17] Lemma 2.2 holds for PH 1 fit, IP) regular vector fields 
when the jump conditions are met on the interfaces. 

Theorem A .6 (Weighted Helmholtz decomposition). Under Assumption A.l, for 
any v E iTo(curl;H), there exist 0 E Hq(Q) and w E PH 1 {LI, HP) n X(fl,/3) such 
that the decomposition (4.2) holds, and the estimate (4.3) is true. 

Proof. For any v E iTo(curl; H), let 0 E Hq{Q) be the solution of 

(0V0, V0) = {/3v, V0), V0 E Hq(LI). 

It is easy to check that 

0 1/2 V0 

and that w = v — V0 satisfies 

V- (fiw) = 0 in Ll and w x n = 0 on dLl. (A.7) 
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L 2 (D.) 


< 




V 


L 2 (Cl 1 


< nil’ll 


(A.6) 










The decomposition v = w + shares the same form of the result (4.2) we want to 
prove, yet the rest is to show that w E PiT 1 (n, &). To this end, we first construct 
an iT 1 -lifting of the w. Using integration by parts we have 

-[ Vcj)- (wxn)dS = [ Vf Vxto = f </>VxwndS V0 E H 1 ^), 

Jon Jn Jan 

thus w x n = 0 on dQ implies V x w ■ n = 0 on Sfl by a density argument (e.g. see 
[1]). Applying Theorem 3.17 in [1] on Vxid, there exists a wo E X (Q, 1) such that 

’ V x wo = V x w in U, 

< V-wo = 0 in U, 

wq x n = 0 on dFl. 

Taking the convexity of the into account, an embedding result from Theorem 2.17 
in [1] reads that X(fl, 1) ^ Thus wq E iT 1 (U). Obviously, 

Vx (w — wo) = 0 in 0 and (w — wq) xn = 0 on dfl. 

The simply-connectedness of O implies that there exists a ( £ H 1 (fi) (see Lemma 
A.2) with a constant boundary value such that 


w — wo = V£ in U. 

By the fact that w | is divergence free within each flj respectively, and \_f3w ■ = 

0 for any F C 3, one can check that the variational problem that Q satisfies is 

(0VC,V0) =-J2 [ Wwo-ni^dS V</> E 
Fc3^ F 

Noticing \J3wq • ri^ F = [/3] F (wo • n)| E H 1 / 2 (F) on any F C 3, the regularity 
result of Theorem 4.1 in [17] shows that £ E PH 2 (fl, ^), in which a function is 
piecewisely H 2 smooth, while has FL 1 regularity across the interfaces on the whole 
domain. This, in turn, implies that 

w = Wo + vc e PH 1 (n 1 ^>)n{v e x(n,p) ■. v-(pv) = o}. 


Lastly, to prove the estimate, by the triangle inequality and (A. 6 ), we have 


/ 3 1 / 2 


w 


L 2 (Q) 


< 


P 1 ' 2 


V 


L 2 (fi) 


+ 


p 1/2 vpj 


L 2 (n) 


< C'|||«|| 


It follows from Assumption A.l (ii) and Lemma A.4 that 


E 

3=1 




-V 2 V 


w 


i 2 (n i ) 


< 


Cj2\\p 1/2 Vw 


3 =1 


L 2 (0. j ) 


= c 


pV 2 V 


x w 


L 2 (p) 


< C'|||«|||. 


These inequalities and (A. 6 ) imply the validity of (4.3) and, hence, it completes the 
proof of the theorem. □ 
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Remark A.7. The decomposition result in Theorem A.6 resembles that of Theorem 
3.5 in [17]: any vector field in X(D,/3) can be split into a PH l (Q, 2?)-regular part, 
and a singular part solving a Dirichlet boundary problem —V-(/3V?/>) = / 6 L 2 (fl). 
In the proof of Theorem A. 6 , we refine the results to cater the need for the pipeline 
of proving the reliability of the error estimator. Namely, when certain assumption of 
geometry is imposed, if a vector field v e X(fi, /?) with its tangential trace vanishing 
on the boundary, and V -(/3v) = 0, that singular part is non-existent. 
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